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On the limit between short and long GRBs 


M. Tarnopolski 


Abstract Two classes of GRBs have been identified 
thus far without doubt and are prescribed to different 
physical scenarios - NS-NS or NS-BH mergers, and col¬ 
lapse of massive stars, for short and long GRBs, respec¬ 
tively. The existence of two distinct populations was 
inferred through a bimodal distribution of the observed 
durations Tgo, and the commonly applied 2 s limit be¬ 
tween short and long GRBs was obtained by fitting a 
parabola between the two peaks in binned data from 
BATSE IB. Herein, by means of a maximum likelihood 
(ML) method a mixture of two Gaussians is fitted to the 
datasets from BATSE, Swift, BeppoSAX, and Fermi in 
search for a local minimum that might serve as a new, 
more proper, limit for the two GRB classes. It is found 
that Swift and BeppoSAX distributions are unimodal, 
hence no local minimum is present, Fermi is consistent 
with the conventional limit, whereas BATSE gives the 
limit significantly longer (equal to 3.38 ± 0.27 s) than 
2 s. These new values change the fractions of short 
and long GRBs in the samples examined, and imply 
that the observed Tgg durations are detector depen¬ 
dent, hence no universal limiting value may be applied 
to all satellites due to their different instrument specifi¬ 
cations. Because of this, and due to the strong overlap 
of the two-Gaussian components, the straightforward 
association of short GRBs to mergers and long ones to 
collapsars is ambiguous. 
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1 Introduction 


Gamma-ray bursts (GRBs) were detected by military 
satellites Vela in late 1960’s. GRBs were recognized 


early to be of extrasolar origin (Klebesadel, Strong & 


Olson 1973). Mazets et al. (1981) first observed a bi¬ 


modal distribution of Tgg (time during which 90% of 
the burst’s fluence is accumulated) drawn for 143 events 
detected in the KONUS experiment. Burst and Tran¬ 
sient Source Explorer (BATSE) onboard the Gomp- 


ton Gamma Ray Observatory (CGRO) (Meegan et al. 


19921 allowed to confirm the hypothesis that GRBs are 


of extragalactic origin due to isotropic angular distri¬ 
bution in the sky combined with the fact that they ex¬ 
hibited an intensity distribution that deviated strongly 


from the —3/2 power law (Briggs 1995 Fishman & Mee¬ 


gan 


19951. However, a more complete sample of BATSE 


short GRBs were shown to be distributed anisotrop- 


ically (Meszaros & Stocek 


2003 


Vavrek et al. 


2008) 


and cosmological consequences were discussed lately 


(Meszaros et al. 2009). BATSE IB data release was 


followed by further investigation of the Tgg distribution 


(Kouveliotou et al. 1993) that lead to establishing the 


common classification of GRBs into short (Tgg < 2 s) 
and long (Tgg > 2 s). This 2 s limit was derived by 
fitting a parabola to the local minimum of the binned 
distribution of 222 GRBs. It was observed that dura¬ 
tions Egg seem to exhibit log-normal distributions which 


were fitted to short and long GRBs (McBreen et al. 


1994), resulting in mean durations equal to 0.37 s and 


26.36 s. A mixture of Gaussians fitted to log Tgg dataset 
from BATSE 2B yielded locations of the components 


equal to 0.60s and 32.1s (Koshut et al. 1996), while a 


subset of BATSE 3B sample yielded 0.42 s and 34.4 s 
( Kouveliotou et ^[1996 ). A complete BATSE dataset 
gave mean locations of the groups at 0.78 s and 34.7 s 
( Horvath||2002 ). The progenitors of long GRBs are as¬ 
sociated with supernovae ( Woosley Sz Bloom]|2006 ) re¬ 
lated with collapse of massive, e.g. Wolf-Rayet, stars. 
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Progenitors of short GRBs are thought to be NS-NS 


or NS-BH mergers (Nakar 20071, and no connection 


between short GRBs and supernovae has been proven 
( Zhang et al.||2009 ). 

The existence of an intermediate-duration GRB 
class, consisting of GRBs with Tqq in the range 2 — 10 s, 
was put forward (Horvath 1998 Mukherjee et aHl998 ) 
based on the analysis of BATSE 3B data. It was sup¬ 


ported (Horvath 20021 with the use of the complete 


BATSE dataset. Evidence for a third normal com- 


2008 


2009 


Zhang & Choi 

2008 

Huja, Meszaros & Ri'pa 

Horvath et al. 2010). 

BeppoSAX dataset was 


shown to be in agreement with earlier results regard¬ 
ing the bimodal distribution, and the detection of an 
intermediate-duration component was established on 
a lower, compared to BATSE and Swift, significance 
level due to a less populate sample ( Horvath]|2009 ). It 
is important to note that in BeppoSAX only the inter¬ 
mediate and long GRBs were detected, the short ones 


being not present. Interestingly, Zitouni et al. (2015) 


re-examined the BATSE current catalog as well as the 
Swift dataset, and found that a mixture of three Gaus- 
sians hts the Swift data better than a two-Gaussian, 
while in the BATSE case statistical tests did not sup¬ 
port the presence of a third component. Regarding 
Fermi, a three-Gaussian is a better fit than a two- 
Gaussiar0 however the presence of a third group in the 
Ego distribution was found to be unlikely (Tarnopolski 
|2015a|bD . 

The 2 s limit is widely used in GRB analysis. How¬ 


ever, the Swift data were re-examined (Bromberg et al. 


2013) and it was found that a limit of 0.8 s is more 


suitable for the GRBs observed by Swift. Many works 
in which a two-Gaussian was fitted to the logTgg dis¬ 
tribution showed a signihcant overlap of components 


corresponding to short and long GRBs (McBreen et 


aL||1994t IKoshut et al. 19961 Horvath 2002 |Zhang & 

odl 

o 

o 

‘o 

O 

Huja, Meszaros & Ri'pa|2009l 

Horvath|2009 


datasets from BATSE, Swift, BeppoSAX, Fermi, among 
others. The mentioned datasets consist of ~ 1000—2000 
events. Based on the well-established conjecture that 
durations Egg are log-normally distributed, the limit 
between short and long GRBs may be placed at the 
position of the local minimum of a mixture distribu¬ 
tion. 


The aim of this paper is to examine what limits are 
most suitable for GRB samples observed by different 
satellites. In Section]^ the datasets and methods are 
described. Results are shown in Section while Sec¬ 
tion]^ is devoted to discussion and is followed by con¬ 
cluding remarks gathered in Section 


2 Datasets and methods 


The dataset^ from BATShj^ BeppoSAi^ and 

Ferm^aie considered herein. They contain 2041 GRBs 
(BATSE current catalog), 914 (Swift), 1003 {Bep¬ 
poSAX), and 1596 (Fermi). Additionally, a subset of 
BATSE data, i.e. complete BATSE IB sample con¬ 
taining 226 GRBs, is examined to compare with results 


of Kouveliotou et al. (1993). Eor display purposes, his¬ 


tograms are plotted using the Knuth rule for bin width. 
Up to date, to the best of the author’s knowledge, only 


Horvath et al. (2012) and Qin et al. (2013) conducted 


research on a Fermi subsample, consisting of 425 GRBs 
from the first release of the catalog. 

The fittings are performed using the maximum like¬ 


lihood (ML) method (Kendall & Stuart 1973). Hav¬ 
ing a distribution with a probability density function 
(PDF) given by / = f(x; 6) (possibly a mixture), where 
9 = is a set of p parameters, the log-likelihood 

function is defined as 


N 


^ = ^^ogf(x^-,e). 


( 1 ) 


where {xi}^i are the datapoints from the sample to 
which a distribution is htted. The fitting is performed 
by searching a set of parameters 9 for which the log- 
likelihood C is maximized. The fitted function in this 
case is a mixture of two standard Gaussians, J\f(fj,,cr^): 


fk(x) ='^A,<p 


X- 


= E 


A,; 


y/^ai 


exp 


[x - 

2a- 

( 2 ) 


Here, = 2, so the distribution is described by p = 5 
parameters: two means /xi, /X 2 , two dispersions fii, <72, 
and one weight Ai. The second weight is ^2 = 1 — 


^Adding parameters to a model always results in a better fit (in 
the sense of a lower or a higher maximum log-likelihood) due 
to more freedom given to the model to follow the data. The 
important question is whether this improvement is statistically 
significant, and whether the model is an appropriate one. See 
( Tarnopolski|2015a|b I for a discussion. 


^All accessed on April 29, 2015. 

"http://gammaray.msf c.nasa.gov/batse/grb/catalog/current 
^ http://swift.gsf c.nasa.gov/archive/grbtable 
"https://heasarc.gsf c.nasa.gov/docs/sax/sax.html 

^ http://heasarc.gsf c.nasa.gov/W3Browse/fermi/fermigbrst. 
html 
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due to normalization. Normal distribution’s PDF is 
denoted by (p. 

To estimate the parameter errors, 86, a simple Monte 
Carlo technique called a parametric bootstrap (|Efron 


1979 1981 Efron k, Tibshirani|1981 1 is performed, i.e., 
having a distribution from Eq. fitted, it is randomly 
sampled to create a set of N random variates {N being 
the same as in the original dataset). This set is used 
to find another fit. After repeating this procedure 1000 
times, the standard deviations are computed from the 
1000 sets of 5 parameters, and serve as errors for the 
parameters obtained from the original dataset. 

In the same manner the error of the location of the 
local minimum is estimated. In this case, some of the 
realisations drawn from the original bimodal distribu¬ 
tion may happen to be unimodal, hence no local mini¬ 
mum might be present. In that case, only a fraction of 
the 1000 realisations which do have a local minimum is 
taken into account. The opposite situation may also oc¬ 
cur, i.e. sampling a unimodal distribution and execut¬ 
ing the bootstrap may result in some bimodal realisa¬ 
tions. However, it turnt out that this kind of situations 
happen rarely, especially the latter, so no ambiguity is 
encountered. 


3 Results 


The results in graphical form are displayed in Fig. 
where the vertical solid line marks the conventional 
limit of 2 s, and the vertical dashed line marks the lo¬ 
cation of the minimum of a mixture of two normal dis¬ 
tributions (a two-Gaussian). The Swift and BeppoSAX 
distributions are unimodal, so no new limit may be in¬ 
ferred. Among the 1000 bootstrap executions, 163 re¬ 
alisations were bimodal for the Swift sample, and only 
one yielded a local minimum in the case of BeppoSAX. 
Parameters of the fits are gathered in Table 

BATSE IB yielded a minimum at 2.16 ± 0.05 s, close 
to the value attained by Kouveliotou et ah] ( |1993 ). 
However, for the current BATSE catalog a limit of 
3.38±0.27 s is more suitable, and the conventional value 
of 2 s lies outside the interval more than five times of 
the error. Fermi dataset is the most consistent with 
the 2 s limit, yielding a minimum at 2.05 ± 0.25 s. The 
shallower the minimum, the bigger the error obtained. 

In the case of BATSE current catalog, all 1000 boot¬ 
strap realisations were bimodal; in BATSE IB a min¬ 
imum was present in 911 cases, while for Fermi there 
were 896 bimodal realisations. 


4 Discussion 


Bromberg et al. (2013) found, by constructing a PDF 


for collapsars and non-collapsars (a classification based 
on physical origin of a GRB), that while the 2 s criterion 
is suitable for BATSE, a value of 0.8 s is more appropri¬ 
ate for the Swift dataset. Herein, the logTgo distribu¬ 
tions examined imply that the suitable limit between 
short and long GRBs for the BATSE current catalog 
should be 3.38 s, based on a univariate analysis. This is 
significantly higher than the commonly applied 2 s cri¬ 
terion. In case of Swift, the duration distribution turnt 
out to be unimodal, and as such no natural limit may 
be inferred. Also BeppoSAX durations are unimodal, 
giving no new limiting value. It is important to note 
that the locations fii, corresponding to the shorter com¬ 
ponent, are negative (hence Tgo < Is) for BATSE (IB 
and current). Swift and Fermi, while for BeppoSAX it is 
pi = 0.626, corresponding to Tgg « 4.23 s. This is defi¬ 
nitely not a short GRB group, and it is consistent with 


(Horvath 2009) where the short GRB group was not 


detected. The Fermi data have a minimum at 2.05 s, 
consistent with the common limiting value. 

The newly obtained limits result in different popula¬ 
tions of short and long GRBs in the datasets examined 
(see Table . In BATSE IB the fraction of long GRBs 
in the sample is unchanged (due to smallness of the 
sample; it appears there are no GRBs with durations 
between 2 s and 2.16 s). In Fermi this fraction is nearly 
the same, slightly smaller than conventional due to the 
limit being slightly higher than 2 s. The biggest differ¬ 
ence is visible in the BATSE current catalog, where the 
new limit leads to diminishing the long GRBs fraction 
by 4%. 

All catalogs are dominated by long GRBs, the high¬ 
est fraction of more than 90% being observed by Swift. 
In BATSE current, the proportion of short and long 
GRBs is ^ 1 : 3. Swift is more sensitive in soft bands 
(corresponding to long GRBs) than BATSE was, while 


Table 2 Fractions of long GRBs and overlap of compo¬ 
nents of the two-Gaussian fits. 


Label 

Dataset 

Long GRBs fraction [%] 

Overlap [%] 

GonventionaF 

This work'’ 

(a) 

BATSE IB 

73.89 

73.89 

5.68 

(b) 

BATSE current 

75.50 

71.53 

10.1 

(c) 

Swift 

90.81 

— 

9.30 

(d) 

BeppoSAX 

88.24 

— 

34.4 

(e) 

Fermi 

83.40 

83.08 

9.06 


“When the conventional | |Kouveliotou et al. 11993] ) limit of 2 s is applied. 
*’When the new limits obtained herein are applied. 
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Table 1 Parameters of the fits. Label corresponds to panels in Fig. Errors are estimated using the bootstrap method. 


Label 

Dataset 

N 

i 


3pi 

(7i 

5(ji 

Ai 

5A, 

min. 

Jmin. 

(a) 

BATSE IB 

226 

1 

-0.393 

0.099 

0.465 

0.069 

0.272 

0.040 

2.158 

0.049 

2 

1.460 

0.056 

0.532 

0.044 

0.728 

(b) 

BATSE 

2041 

1 

-0.095 

0.051 

0.627 

0.033 

0.336 

0.018 

3.378 

0.272 

current 

2 

1.544 

0.018 

0.429 

0.013 

0.664 

(c) 

Swift 

914 

1 

-0.026 

0.255 

0.740 

0.120 

0.139 

0.042 



2 

1.638 

0.031 

0.528 

0.023 

0.861 



(d) 

BeppoSAX 

1003 

1 

2 

0.626 

1.449 

0.186 

0.035 

0.669 

0.393 

0.075 

0.027 

0.355 

0.645 

0.084 

— 

— 

(e) 

Fermi 

1596 

1 

-0.072 

0.073 

0.525 

0.044 

0.215 

0.021 

2.049 

0.248 

2 

1.451 

0.021 

0.463 

0.014 

0.785 



logTgo 


Fig. 1 Two-Gaussian PDFs fitted to logTgo data. Color dashed curves are the components of the (solid black) mixture 
distribution. Vertical solid line marks the conventional 2 s limit between short and long GRBs. Vertical dashed line marks 
the position of the local minimum (if present) of the mixture. The panels correspond to (a) BATSE IB, (b) BATSE current, 
(c) Swift, (d) BeppoSAX, and (e) Fermi catalogs. In the latter, the new limit is very close to the conventional limit of 
Kouveliotou et ^ (19931. 
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Fermfs sensitivity at very soft and very hard GRBs had 19931, whereas in BATSE current catalog the value ob- 


increased compared to BATSE (Meegan et al. 20091. 


BeppoSAX is also more sensitive to long GRBs due to 
the trigger system which used 1 s as short integration 
time ( Horvath||2009| ) , hence the lack of a distinct short 
GRB peak. GRBs tend to be softer at later times, hence 
the inferred duration is shorter than it might be. This 
naturally leads to a conclusion that the duration dis¬ 
tributions as observed by different satellites must differ 
between each other, and also the limit between short 
and long GRBs (the local minimum) has to be placed 
at different locations. 

The duration Tgg itself is not an unambiguous indi¬ 
cator of a GRB type, as the components of the fitted 
two-Gaussians overlap strongly. To quantify this over¬ 
lap, the common area under the curves is computed (the 
total area of a two-Gaussian PDF is equal to unity, and 
the area under each component is given by the weights 
Ai). This gives a probability of misclassifying a GRB 
from 5.68% (BATSE IB) to 10.1% (BATSE current), 
and an enormous 34.4% in the case of BeppoSAX (see 
Table [§. 

A solution, proposed to deal with the classification 
ambiguity problem, was proposed and examined in a 


number of papers (Hakkila et al. 

2003 Horvath et al. 

2004 

2006[ Chattopadhyay et a 

.| 20071 Veres et al. 

2010 

Horvath et al. 2010). The idea is to examine 


a multi-dimensional space of various parameters; par¬ 
ticularly, a two-dimensional space of the hardness ra¬ 
tio vs. duration Tgg. This approach still awaits to 
be applied to the Fermi GRBs. Additional parame¬ 
ters have been defined and proposed for GRB classifi¬ 
cation as well. Examples are e = ,iso,52/E’^.'1^2 (™- 


ambiguously dividing short and long GRBs) (Lii et al. 


2010 

), minimum variability time-scale (MVTS) (Bhat 

2013 

MacLachlan et al.|2012[ 

2013a|bl|Golkhou & But- 

ler||2014 Golkhou, Butler & Littlejohns|2015 ) or Hurst 


2015c I. Still, the most common criterion is the GRB 


duration, and its limitting value has been shown herein 
to be detector dependent. 


tained is equal to 3.38 ±0.27 s. This leads to a different 
than commonly established, fraction of long GRBs in 
the sample, diminished by 4% (see Table [^. 

Due to the significant overlap and dependence of the 
location of the minimum on the detector, while the divi¬ 
sion into short and long GRBs based on their durations 
is qualitatively proper, it is not unambiguously related 
to its progenitor, i.e. collapsar or non-collapsar. There¬ 
fore, as the short-long phenomenological classification 
justifies the existence of two distinct GRB classes, it 
gives limited insight into the underlying physical phe¬ 
nomenon. 


5 Conclusions 

The duration distributions of various catalogs (BATSE 
IB, BATSE current. Swift, BeppoSAX, and Fermi) were 
examined. A mixture of two Gaussians was fitted to 
the logTgo distributions in search for a new limiting 
value placed at the local minimum. It was found that 
the datasets from Swift and BeppoSAX are unimodal, 
hence no new limit may be inferred. The results from 
BATSE IB and Fermi are consistent with the conven¬ 


tional phenomenological limit of 2 s (Kouveliotou et al. 



























































References 

Barnacka A., Loeb A.: Astrophys. J.794, L8 (2014) 

Bhat P. N.: in Castro-Tirado A. J., Gorosabel J., Park I. 
H., eds, EAS Publ. Ser. Vol. 61, Temporal Decomposi¬ 
tion Studies of GRB Lightcurves. Cambridge Univ. Press, 
Cambridge, p. 45 (2013) 

Briggs M. S.: Astrophys. Space Sci.231, 3 (1995) 

Bromberg O., Nakar E., Piran T., Sari R.: Astrophys. J.764, 
179 (2013) 

Chattopadhyay T., Misra R., Chattopadhyay A. K., Naskar 
M.: Astrophys. J.667, 1017 (2007) 

Efron B.: Ann. Stat. 7, 1 (1979) 

Efron B.: Biometrika 68, 589 (1981) 

Efron B., Tibshirani R. J.: An Introduction to the Boot¬ 
strap. CRC Press (1994) 

Fishman G. J., Meegan C. A.: Annu. Rev. Astron. Astro- 
phys.33, 415 (1995) 

Golkhou V. Z., Butler N. R.: Astrophys. J.787, 90 (2014) 
Golkhou V. Z., Butler N. R., Littlejohns O. M.: Astro¬ 
phys. J.in press (arXiv:1501.05948) (2015) 

Hakkila J., Giblin T. W., Roiger R. J., Haglin D. J., Paciesas 
W. S., Meegan C. A.: Astrophys. J., 582, 320 (2003) 
Horvath L: Astrophys. J.508, 757 (1998) 

Horvath L: Astron. Astrophys.392, 791 (2002) 

Horvath L, Meszaros A., Balazs L. G., Bagoly Z.: Balt. 
Astron., 13, 217 (2004) 

Horvath I., Balazs L. G., Bagoly Z., Ryde F., Meszaros A.: 

Astron. Astrophys.447, 23 (2006) 

Horvath L, Balazs L. G., Bagoly Z., Veres P.: Astron. As¬ 
trophys.489, LI (2008) 

Horvath L: Astrophys. Space Sci.323, 83 (2009) 

Horvath L, Bagoly Z., Balazs L. G., de Ugarte Postigo A., 
Veres P., Meszaros A.: Astrophys. J.713, 552 (2010) 
Horvath L, Balazs L. G., Hakkila J., Bagoly Z., Preece R. 

D.: PoS (GRB2012), 046 (2012) 

Huja D., Meszaros A., Ripa J.: Astron. Astrophys.504, 67 
(2009) 

Kendall M., Stuart A.: The Advanced Theory of Statistics. 
Griflin, London (1973) 

Klebesadel R. W., Strong 1. B., Olson R. A.: Astro¬ 
phys. J.182, L85 (1973) 

Koshut T. M., Paciesas W. S., Kouveliotou C., van Paradijs 

J. , Pendleton G. N., Fishman G. J., Meegan G. A.: As¬ 
trophys. J.463, 570 (1996) 

Kouveliotou C., Meegan G. A., Fishman G. J., Bhat N. P., 
Briggs M. S., Koshut T. M., Paciesas W. S., Pendleton 
G. N.: Astrophys. J.413, LlOl (1993) 

Kouveliotou C., Koshut T., Briggs M. S., Pendleton G. N., 
Meegan C. A., Fishman G. J., Lestrade J. P.: AIP Conf. 
Proc. 384, 42 (1996) 

Lii H.-J., Liang E.-W., Zhang B.-B., Zhang B.: Astro¬ 
phys. J.725, 1965 (2010) 

MacLachlan G. A., Shenoy A., Sonbas E., Dhuga K. S., 
Eskandarian A., Maximon L. C., Parke W. G.: Mon. Not. 
R. Astron. Soc.425, L32 (2012) 

MacLachlan G. A. et ah: Mon. Not. R. Astron. Soc.432, 857 
(2013a) 

MacLachlan G. A., Shenoy A., Sonbas E., Coyne, R., Dhuga 

K. S., Eskandarian A., Maximon L. C., Parke W. C.: 
Mon. Not. R. Astron. Soc.436, 2907 (2013b) 


Mazets E. P., Golenetskii S. V., Ilinskii V. N., Panov V. N., 
Aptekar R. L., Gurian I. A., Proskura M. P., Sokolov 1. 

A. , Sokolova Z. L, Kharitonova T. V.: Astrophys. Space 
Sci.80, 3 (1981) 

McBreen B., Hurley K. J., Long R., Metcalfe L.: Mon. Not. 

R. Astron. Soc.271, 662 (1994) 

Meegan C. A., Fishman G. J., Wilson R. B., Horack J. M., 
Brock M. N., Paciesas W. S., Pendleton G. N., Kouve¬ 
liotou G.: Nature355, 143 (1992) 

Meegan C. A., Lichti G., Bhat P. N., Bissaldi E., Briggs 
M. S., Gonnaughton V., Diehl R., Fishman G., Greiner 
J., Hoover A. S., van der Horst A. J., von Kienlin A., 
Kippen R. M., Kouveliotou C., McBreen S., Paciesas W. 

S. , Preece R., Steinle H., Wallace M. S., Wilson R. B., 
Wilson-Hodge, C.: Astrophys. J.702, 791 (2009) 

Meszaros P., Stocek J.: Astron. Astrophys.18, 293 (2003) 
Meszaros A., Balazs L. G., Bagoly Z., Veres P.: Balt. As¬ 
tron. 18, 293 (2009) 

Mukherjee S., Feigelson E. D., Jogesh Babu G., Murtagh 
F., Fraley C., Raftery A.: Astrophys. J.508, 314 (1998) 

Nakar E.: Phys. Rep.442, 166 (2007) 

Qin Y., Liang E.-W., Liang J.-F., Yi S.-X., Lin L., Zhang 

B. -B., Zhang J., Lii H.-J., Lu R.-J., Lii L.-Z., Zhang B.: 
Astrophys. J.763, 15 (2013) 

Tarnopolski M.: Astron. Astrophys.in press (arXiv: 1506.07324)| 
(2015a) 

Tarnopolski M.: preprint (arXiv: 1506.07801) (2015b) 
Tarnopolski M.: preprint (arXiv: 1507.04886) (2015c) 

Woosley S. E., Bloom J. S.: Annu. Rev. Astron. Astro- 
phys.44, 507 (2006) 

Vavrek R., Balazs L. G., Meszaros A., Horvath L, Bagoly 
Z.: Mon. Not. R. Astron. Soc., 391, 1741 (2008) 

Veres P., Bagoly Z., Horvath L, Meszaros A., Balazs L. G.: 

Astrophys. J.725, 1955 (2010) 

Zhang Z.-B., Choi C.-S.: Astron. Astrophys.484, 293 (2008) 
Zhang B., Zhang B.-B., Virgili F. J., Liang E.-W., Kann 
D. A., Wu X.-F., Proga D., Lv H.-J., Toma K., Meszaros 
P., Burrows D. N., Roming P. W. A., Gehrels N.: Astro¬ 
phys. J.703, 1696 (2009) 

Zitouni H., Guessoum N., Azzam W. J., Mochkovitch R.: 
Astrophys. Space Sci.357:7 (2015) 


This manuscript was prepared with the AAS lATj^X macros v5.2. 




